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We have used soft x-ray magnetic diffraction at the Fe 3+ L2,3 edges to examine to what extent 
the Dzyaloshinsky-Moriya interaction in Ba3NbFe3Si2 0i4 influences its low temperature magnetic 
structure. A modulated component of the moments along the c-axis is present, adding to the 
previously proposed helical magnetic configuration of co-planar moments in the a, 6-plane. This 
leads to a "helical- butterfly" structure and suggests that both the multi-axial in- plane and the 
uniform out-of- plane Dzyaloshinsky-Moriya vectors are relevant. A non zero orbital magnetic signal 
is also observed at the oxygen K edge, which reflects the surprisingly strong hybridization between 
iron 3d and oxygen 2p states, given the nominal spherical symmetry of the Fe 3+ half filled shell. 



I. INTRODUCTION . 

The term chirality was first utilized in science by Lord ! 
Kelvin. His original definition has evolved with time and ' 
we now speak about a chiral system if such a system ' 
exists in two distinct (enantiomeric) states that are in- ' 
terconverted by space inversion, but not by time rever- f 
sal combined with any proper spatial rotation. 1 Chiral- f 
ity permeates natural sciences from biochemistry to solid ' 
state physics. The fact that living organisms use only the f 
left enantiomers of amino acids is still not well under- ' 
stood. Chirality is also found in magnets. 2,3 An example f 
is the left- or right- handedness associated with the he- ' 
lical order of magnetic moments. In principle, the two f 
states are degenerate, resulting in an equipopulation of f 
chiral domains. However, competing interactions or ex- ' 
ternal effects such as strain, can unbalance this ratio, ' 
favoring one particular state. In particular, in non cen- ' 
trosymmetric crystals, characterized by the absence of' 
parity symmetry, a single domain might be selected. De- ' 
spite having 65 non centrosymmetric (including 22 chi- ' 
ral) space groups allowing chiral crystal structures, out ' 
of 230, only few single handed magnetic compounds were ' 
reported. 4-7 Interest in such systems is two-fold. First, ' 
they can exhibit interesting physical properties such as ' 
magnetic Skyrmion lattices 8 or helimagnons. 7 The sec- ' 
ond is related to the discovery of magnetically induced ! 
multiferroics 9 where researchers struggle to find mate- ' 
rials with a stronger electrical polarization. 10 The lat- ! 
ter is directly affected by the imbalance between chiral ' 
domains, which possess opposite electric polarizations. ' 
Therefore, materials showing a single chiral domain are 1 
promising candidates to host a significant macroscopic 1 
electrical polarization, which makes them an ideal model 1 
system to study. Ba3NbFe3Si20i4 gathered attention in ! 
this respect, exhibiting fully chiral magnetism 5 and mag- 1 
nctoclcctric coupling phenomena. 



Ba3NbFc3Si2 0i4 crystallizes in a trigonal P321 space 
group (a = b = 8.539, c = 5.241, 7 = 120°). It displays 
an antifcrromagnctic order below Tjv=27 K. The mag- 
netic moments are localized on the Fe 3+ ions (L ~ 0, 
S = 5/2). These occupy the Wyckoff position (3f) 
(0.2496, 0, 0.5) with .2. site symmetry, forming trian- 
gular units in the a,6-planes. Elastic neutron scattering 
studies 5 suggest that the same triangular configuration 
of co-planar moments at 120° from each other is stabi- 
lized within each triangle of an a, 6-plane and that this 
arrangement is helically modulated from a, 6-plane to a, 6- 
plane along the c-axis according to the propagation vec- 
tor (0, 0, t) with r close to 1/7 (see Fig. 1 of Ref. 5). An 
extremely appealing discovery was that the single crys- 
tals are grown enantiopure and that the low tempera- 
ture magnetic structure is single domain, with a single 
chirality of the triangular magnetic arrangement on the 
triangles and a single chirality of the helical modulation 
of the magnetic moments, which was dubbed helicity. 5 
It was suggested that the Dzyaloshinsky-Moriya 14,15 ex- 
change interaction might be responsible for selecting the 
ground state configuration 5 and for the opening of a small 
gap in the magnetic excitation spectrum. 16 Another in- 
elastic neutron scattering study proposed the latter to 
arise from single ion anisotropy, 17 but recent spin reso- 
nance experiments support the first scenario indicating 
furthermore that not only the uniform component along 
the c-axis of the Dzyaloshinsky-Moriya vector but also 
its multiaxial component within the a, 6-plane might be 
sizeable. 18 The latter could generate an additional com- 
ponent to the magnetic structure not necessarily detected 
by neutron scattering. To find evidence for such a mag- 
netic motif wc have used resonant x-ray diffraction at the 
Fe L edges. Our results show clear deviations from the 
magnetic structure previously proposed, confirming the 
existence of such a component. 
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so II. EXPERIMENTAL DETAILS i« 

145 

91 Powders of Ba 3 NbFe3Si20i 4 were synthesized by solid 146 

92 state reaction from stoichiometric amounts of Nb2C>3, 147 

93 FC2O3, SiC>2 oxides and BaCC>3 barium carbonate, at 148 

94 1150° C in air within an alumina crucible. The reagents 149 

95 were carefully mixed and pressed at lGPa to form com- 

96 pact cylinders before annealing. The phase purity was 

97 checked by x-ray powder diffraction. Single crystals were 

98 grown from the as-prepared polycrystalline cylinders by 150 

99 the floating-zone method in an image furnace. 19 The sin- 151 

100 gle crystal used in the present investigation was extracted 152 

101 from the same batch as the one used in Ref. 5 and has the 153 

102 same structural chirality ct, to be precise ct = — 1- After 154 

103 polishing the surface perpendicular to the [001] direction 155 

104 it was annealed to improve the surface quality. 156 

105 We have performed resonant x-ray diffraction exper- 157 

106 iments at the Fc L2.3 edge. These energies correspond 158 

107 to a wavelength of approximately 17 A and are asso- 159 
10s ciated to an electric dipole resonance from the iron 2p 160 
109 to 3d levels. Experiments were performed with the RE- 161 
no SOXS chamber 20 at the XI IMA beamline 21 of the Swiss 162 
in Light Source. The twin Apple undulators provide linear, 163 

112 horizontal tt and vertical a, and circularly, right R and 164 

113 left L, polarized x rays with a polarization rate close to 165 

114 100%. The polarization of the diffracted beam was not 166 

115 analyzed. The sample was attached to the cold finger 167 
lie of an He flow cryostat with a base temperature of 10 K. 168 
117 Azimuthal scans were achieved by rotation of the single 169 
us crystal, with an accuracy of approximately ±5°. 170 



119 III. RESONANT X-RAY SCATTERING 173 

174 

120 The x-ray cross section for magnetic scattering is nor- 175 

121 mally very small, though at synchrotron photon sources 

122 such weak signals are routinely measurable. However, 

123 when working close to an atomic absorption edge the 

124 magnetic scattering signals are significantly enhanced 

125 and are element sensitive. Resonant x-ray diffraction 178 

126 occurs when a photon excites a core electron to empty 

127 states, and is subsequently re-emitted when the electron™ 

128 and the core hole recombinc. 22-24 This process introducesiso 

129 anisotropic contributions to the x-ray susceptibility ten-isi 



130 sor 
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the amplitude of which increases dramatically™ 
as the photon energy is tuned to an atomic absorption™ 
edge. In the presence of long-range magnetic order, or™ 
a spatially anisotropic electronic distribution, the inter-™ 
ference of the anomalous scattering amplitudes may lead™ 
to Bragg peaks at positions forbidden by the crystallo-™ 
graphic space group. An example of such a resonant™ 
enhancement of the diffracted intensity as a function of™ 
energy occurring in the vicinity of the Fe L3 edge in™ 
Ba 3 NbFe 3 Si20i4 is given in Fig. 1. X rays thus provem 
to be a valid alternative or complementary tool to neu-wa 
tron diffraction for the study of magnetic structures. 28 ~ 32 ™ 
Its superior resolution in reciprocal space can be advanta-™ 
geous, simplifying for instance the precise determination^? 



of incommensurate magnetic phases, which is relevant in 
cases where the incommensurability is very small. 33 

To understand the content of the x-ray resonant mag- 
netic cross section, it is customary to use the expres- 
sion first derived by Hannon and Trammell for an electric 
dipole (El) event: 22 " 24 

F$% = (e' • e)F<°> - i(e' x e) • i n F^ + (e' • z„) (e • z n )F^ , 

(1) 

where the first term contributes to the charge (Thomp- 
son) Bragg peak. The second and third terms correspond 
to magnetic diffraction. z„ is a unit vector in the direc- 
tion of the magnetic moment of the nth ion in the unit 
cell and e (e') describes the polarization state of the in- 
coming (outgoing) x rays. depend on atomic prop- 
erties and determine the strength of the resonance. 2,24 In 
an antiferromagnet, the second term produces the first- 
harmonic magnetic satellites and the third term, which 
contains two powers of the magnetic moment, produces 
the second-harmonic magnetic satellites. It shows how 
the intensity of the magnetic diffraction depends on the 
motif of the magnetic moments and on the orientation 
of the sample relative to the incident x-ray polarization 
state. In particular, a non collinear magnetic motif is able 
to produce a different diffraction intensity depending on 
the helicity of the incident x rays, e.g. Ir 7^ II, where 
Ir is the intensity measured with incident right-handed 
circularly polarized photons and II for left-handed ones. 
Rotating the sample about the diffraction wave vector 
might result in a smooth change of the diffracted intensity 
which helps to reconstruct the magnetic moment motif. 

It is worth emphasizing that Eq. (1) is an approxi- 
mation for the resonant magnetic scattering cross sec- 
tion which, strictly speaking, is only valid for a cylindri- 
cal symmetrical environment of the resonant ion. When 
this approximation does not hold the diffracted intensi- 
ties must be described as exemplified in Rcf. 2, 34-38. 



IV. RESULTS 

Once the sample is cooled below the Neel temperature 
Tjv, superstructure peaks (0, 0, tit) of order n up to three 
arise from magnetic ordering and magnetically induced 
lattice distortions. The observation of such reflections is 
remarkable as, given the magnetic motif suggested by 
neutron diffraction, they should be absent. They are 
of resonant nature and they disappear when the energy 
of the incident x ray is detuned from the iron L edges 
(Fig. 1). Non-resonant magnetic intensity could be zero 
or too small to be visible. Resonant x-ray diffraction is 
sensitive to the spin, orbital and charge degrees of free- 
dom. 24,39-41 In order to assert their origin and refine the 
magnetic structure, we collected their energy, azimuthal 
and temperature dependence. 

Fig. 1 shows the energy dependence of the three super- 
structural peaks collected for x rays with polarization in 
the diffraction plane (so-called tt geometry) and perpen- 
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FIG. 1. (Color online) Intensity versus energy of the three 
satellite reflections in the vicinity of the Fe L3 edge. Spectra 
collected with incident n [(blue) square] and a polarizations 
[(red) filled circle] at 10 K. Spectra are scaled [(0, 0, r) and (0, 
0, 3r) were multiplied by 2.5 and 80 respectively] and shifted 
for clarity and lines are guides to the eye. The black contin- 
uous line represents the sample absorption spectra collected 
in fluorescent mode. 



FIG. 2. (Color online) Azimuthal angle dependence of the 
Itr/Iir (red) and Il/Ir (black) ratio for the (0, 0, r) magnetic 
reflection. The (black) line represents the predictions of the 
model described in the text (\ 2 = 4.0 for comparison with 
both dataset, \ 2 = 1-5 f° r the ratio Il/Ir alone). Measure- 
ments were performed in the vicinity of the Fe L3 edge. The 
azimuthal angle equals zero when the [f 00] direction is in the 
scattering plane. 



us dicular to it (a geometry). They measure the maximum 2 29 

199 intensity of the diffraction peak at different energies (i.e. 230 

200 energy scans at fixed momentum transfer). The first har-231 

201 monic peak (n = 1) shows equal intensity (1^ = I a ) for 232 

202 both incident x-rays polarization as the energy of the233 

203 incident x-rays is swept across the iron L3 edge. The ra-234 

204 tio I n over I a is very close to one and has no significant235 

205 modulation as the sample is rotated about the diffraction 2 36 

206 wave vector (0, 0, r) (so-called azimuthal-angle rotation), 237 

207 as exemplified in Fig. 2. Data are collected for a Bragg 2 38 

208 angle Or = 14.1° where a significant contribution froni239 

209 specular reflectivity is present. Such a contribution is24o 

210 different for 1^ and I a and, combined with the weakness24i 

211 of the signal, complicates the determination of the mag-242 

212 netic Bragg diffraction contribution. In this respect, the243 

213 data gathered with incident circularly polarized photons244 

214 {Ir and II) provide a more reliable data set, as being a245 

215 complex combination of the linearly polarized light, they 246 

216 present the same background for Ir and II- Indeed the247 

217 ratio II over Ir is very close to one over the investigated248 

218 range and sports smaller error bars. 249 

219 The second harmonic (0, 0, 2t) energy dependence has 250 

220 I v I a . Being associated with small lattice or electron 251 

221 density deformations induced by the magnetic ordering, 252 

222 it is expected to exhibit a la/In ratio different from one. 253 

223 We do not observe any intensity far from the absorp- 254 

224 tion edge. It indicates that the signal originates from the 255 

225 asymmetry of the electron density that appears below the 256 

226 magnetic ordering temperature, possibly triggered by the 257 

227 antiferromagnetic ordering. We have also collected its az-258 

228 imuthal angle dependence (Fig. 3). In analogy with the259 



first harmonic peak it shows no modulation, with I a and 
I v constant within the error bars. Such results are sup- 
ported also by the azimuthal variation of the ratio I a over 
I n which displays smaller error bars due to the elimina- 
tion of possible systematic errors, which equally affect 
both intensities, such as misalignments and changes in 
the sample illuminated area during the azimuthal scan. 
Finally we discuss the third harmonic reflection (0, 0, 
3t). Its energy dependence is quite peculiar. Being /„. 
equal to I a suggests the peak to be of magnetic origin, as 
in the case of (0, 0, r) reflection. However, the spectral 
shape differs strongly from the one of the fundamental 
harmonic. It presents two principal features close in en- 
ergy rather than a single peak with two shoulders as in 
the case of the (0, 0, r). As the iron site symmetry 
(.2.) docs not forbid mixed events (e.g. electric dipolc- 
quadrupole) one possible explanation can be a small con- 
tribution coming from the electric quadrupole or electric 
dipole-quadrupole event, though such contributions are 
usually expected to be negligible. Note that the odd 
reflection intensities are indeed very small compared to 
other magnetic ordering signal found in oxides. 28,30,42-46 
Effect of absorption correction can be discarded as they 
would influence more significantly the (0, 0, t) reflection. 
At lower angles the penetration length is reduced as the 
x rays have to travel longer into the sample before be- 
ing diffracted into the detector. It was unfortunately not 
possible to collect its azimuthal angle dependence due to 
the weakness of the signal. 

The temperature dependence of the satellite reflections 
(Fig. 4) shows strong resemblance to the one observed in 
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FIG. 3. (Color online) Azimuthal angle dependence of the (0, 
0, It) superstructural reflection. The line represents a fit to 
the data with a constant {\ 2 = 1.6 for the ratio 7 (T /7, r ), as 
expected form the model presented in the text. Measurements 
were performed in the vicinity of the Fe L3 edge (E=709.4 eV). 
The azimuthal angle equals zero when the [100] direction is 
in the scattering plane. 



rare-earth metals. 47 Pursuing the parallel with the rare- 
earth metals we would expect that the first harmonic 
arises from magnetic diffraction at the dipole resonance. 
The second harmonic corresponds to charge or orbital 295 
diffraction arising from lattice or electron density modu- 296 
lations. The third-order harmonic might be a magnetic 297 
harmonic of the first or might originate from an elec- 298 
trie quadrupolc resonance, although such a contribution 299 
is expected to be orders of magnitudes weaker. Our es- 300 
timate of the critical exponent (3 found that it is not 301 
consistent with mean-field theory. A fit to power-law 302 
behavior I nT oc (T/v — T) 2 ^ n gave an estimate for the 303 
critical exponents. They are respectively f}\ =0.34±0.04, 304 
j3 2 =0.54±0.05, A3 =0.93±0.08. 

Given the long modulation period of the magnetic 306 

307 

structure it was possible to extend our investigation also 
to the oxygen K edge, which corresponds to an electric 
dipolar transition from the Is to the 2p level. Upon 
cooling below Tjy a signal is observed at this energy. 308 
Figure 5 shows its resonant nature. Observation of a 
resonant signal on an anion is not unusual. 48-50 A reso-309 
nant signal can arise, given a non zero overlap betweemio 
the initial and the final state, whereas a difference existS3ii 
in the up/down spin dipolar overlap integrals. The dif-312 
ference can be induced by polarization of the orbitals. 51 3i3 
Such an asymmetry can arise also in case of a difference3i4 
in the lifetime of the up/down spin channels. Recently3i5 
Beale et al. 49 observed a resonant signal at the oxygen K316 
edge in TbMi^Os, which they interpreted as a signature3i7 
of an antifcrromagnctically ordered spin polarization on3i8 
the oxygen site. Such an observation is quite rcmarkable3i9 
and we share their opinion that the study of oxygen spin32o 
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FIG. 4. Normalized integrated intensity vs temperature of 
the three satellite reflections. The solid lines show the best fit 
to power-law behavior I„ T oc (Tjv - T) 2/3 ". The dashed line 
is the expected mean-field theory dependence. The (0, 0, 2r) 
satellite is 7 times more intense than the (0, 0, r). The same 
ratio holds between the (0, 0, r) and the (0, 0, 3r) satellites. 
The inset (left to right) shows different types of magnetic 
ordering: a simple spiral, a ferromagnetic (conical) spiral, a 
complex spiral (or butterfly) and a static longitudinal wave. 
Data was measured with 7r incident photon energy of 710 eV. 



polarization may lead to new insight in the understanding 
of the magnetoelectric coupling mechanism. As a mat- 
ter of fact, an antiferromagnetic order at the oxygen site 
is consistent with neutron diffraction experiments that 
have already suggested a spin polarization of the oxy- 
gen by finding a value of 4 /is instead of the expected 
5 [Ib for the spherical Fc 3+ half filled ion magnetic mo- 
ment. 5,11 In our case the signal at the oxygen K edge is 
90 times weaker than the corresponding one observed at 
the iron L3 edge. Note that at the K edge the signal 
originates solely from the orbital magnetic moment com- 
ponent. No intensity was observed at the (0, 0, 2r) and 
(0, 0, 3t) satellites at the oxygen K edge. 



V. DISCUSSION 

Insights into the results can be obtained from group 
representation analysis, 52 provided that a single irre- 
ducible representation is selected at the magnetic order- 
ing. The analysis is simplified by the fact that the space 
group P321 associated with the paramagnetic phase is 
symmorphic. It is, to be precise, a semi-direct product 
of the abelian translation group associated with a hexag- 
onal lattice and the dihedral point group 32, which con- 
sists of the identity 1, the anti-clockwise rotation 3 + and 
the clockwise rotation 3~ = (3 + ) 2 about the ternary c- 
axis and the dyads (7r-rotations) about the three binary 
axes at 120° to each other within the a,£>-planc. A vector 
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FIG. 5. (Color online) Intensity [(red) circle] vs energy of the 370 
(0, 0, r) satellite reflection at the oxygen K edge collected at 
10 K with 7T incident x rays. The fluorescence spectra [(blue) 
open square] obtained in the vicinity is also shown. Full (red) 
circle results from a fit of the integrated intensity of a recipro- 373 
cal lattice scan along the c* reciprocal lattice direction. Open 374 
(red) circle are a result of an energy scan with fix momentum375 
transfer. The (black) continuous line is a Gaussian fit of the376 
oxygen resonance with a FWHM = 1.4±0.1 eV. 377 



3 21 along the reciprocal c*-axis is reversed under the dyads 380 

322 and is invariant otherwise. It follows that the star of the 381 

323 magnetic propagation vector consists of the two vectors 382 

324 t± = (0, 0, ±r) each being associated with the little space 383 

325 group P3, which is a semi-direct product of the transla- 384 

326 tion group of the paramagnetic phase and the abelian 385 

327 cyclic point group 3. An abelian group G of no elements 386 

328 has uq conjugacy classes (each being reduced to a single- 387 

329 ton owing to the commutativity) , which implies that it 388 

330 has no irreducible representations Ti (i = 1, ... , tig). It 389 

331 follows that these are necessarily all of dimension d% = l, 390 

332 to comply with the identity Y^i=\ dj = n G- Each Ti co-mi 

333 incides then with its character Xi ■ The value of Xi on any 392 

334 group element g is an nG-th root c 427rp// " G (jp = 1, ... , «g) 393 

335 of 1 , because the order of g always divide . The charac- 394 

336 ter table is then built by making use of the orthogonality 395 

337 theorems. The basis vector of the invariant subspace of 

3 38 each Ti is also easily deduced by applying the projection 

339 operator Vi = ^ J2 g eG Xi(d)*9 011 trial vectors. Table I 
34ii summarizes such results for the cyclic group 3. 

342 The choice of a propagation vector amounts to select 

343 an irreducible representation of the translation group 

344 and determines a dephasing of moments within each 

345 Bravais lattice. Information on the phase relations be- 

346 tween moments of distinct Bravais lattices can be ex- 

347 tractcd only from the irreducible representations of the 

348 little co-group. Three Bravais lattices C u {y = 1, 2, 

349 3) are associated with the positions (0.2496, 0, 0.5), 

350 (0, 0.2496, 0.5) and (-0.2496, -0.2496, 0.5) of the Fe 3 + 

351 ions on the 3f site. Under the symmetry operation 3 + 

352 a moment of C\ (rcsp. £2, £3) is rotated by an an- 



gle of 120° about the c-axis and is transported into £2 
(resp. £3, £1) whereas under the symmetry 3~ it is ro- 
tated by an angle of 240° about the c-axis and is trans- 
ported into £3 (resp. £1, £2). This defines a repre- 
sentation r of the cyclic group 3 of dimension 9 whose 
character \ takes the values x(l) = 9, x(3 + ) = and 
x(3 - ) = on the group elements. T reduces into irre- 
ducible components as : T — 3Fi©3r2©3r3. A magnetic 
structure can be most generally regarded as composed 
of several sine- wave amplitude modulations of moments: 
\{v v {9 v , cv)e~' i ^e~ rr± ' 7V ™ + c.c), where f vn — r v + R n 
defines the position of the moment of C u in the n-th unit 
cell, £„ stands for an initial phase and c.c. means to 
take the complex conjugate. The reduction of T then 
suggests that, whatever the selected irreducible repre- 
sentation Ti, three independent directions of the mo- 
ments are allowed by symmetry and can be combined, 
for instance along two orthogonal unit vectors in the a,b- 
plane, x v = (/t/2, at an angle 0„ from the a-axis and 
Vv = (tt/2, 0jy+7r/2) at an angle 0„ + 7r/2 from the a-axis, 
and along the unit vector z v = (0, 0) of the c-axis, with 
possibly vectors v v {9 Vl cj) v ) of different lengths. 

It was shown, from collected neutron diffraction inten- 
sities, that a hclicoidal modulation is stabilized within 
each £ t/ , associated with a combination of the form 
tv(7r/2,ftv)e- 4 «'- + ae H v u (^,^ + Ti/2)e^^^ = 
rn a .b(&v + iaeHVv)eT % ^" with a = +1 for r+ and a 

I for r_. It is implicitly assumed that the vectors 
v u (ir/2,(/) u ) = m a ^x v and v v {-k j2^ v + n/2) = m a ^y v 
have the same length m a ^, which leads to a circular 
helix. An elliptic helix would have been obtained oth- 
erwise, which a priori cannot be excluded, en = ±1 
defines the magnetic helicity, that is to say the sense of 
the rotation of the moments in the helix as one moves 
along the propagation vector: m(f un ) x fh(f vn + cj = 
e HTia b sin(27rr)(c/| c |) whatever the chosen description 
between r+ and r_. If we impose 0„=2.3 — 4>i according 
to Table I then we must have £1 = £2 = £3, which can 
be set to 0, together with </>i, without loss of general- 
ity. Table I illustrates that a triangular configuration of 
the moments on each triangle is associated with Ti with 
a magnetic triangular chirality +1, that is to say with 
an anti-clockwise sense of the rotation of the moments 
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TABLE I. Irreducible representations of the cyclic point group 
3, little co-group of the propagation vectors f± = (0, 0, ±r) in 
the space group P321, and associated invariant basis vectors. 
v v (6, <f>) symbolizes a vector associated to a Bravais lattice C v 
at an angle 9 from the c-axis and the projection of which in 
the perpendicular plane is at an angle 4> from the a-axis. 
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396 as one moves anti-clockwise on a triangle. A triangu-454 

397 lar configuration of the moments on each triangle witli455 

398 the opposite magnetic triangular chirality —1, that is t0456 

399 say with a clockwise sense of the rotation of the mo- 
no ments as one moves anti-clockwise on a triangle, emerges 

401 from T2 (resp. 1^) when sr = +1 (rcsp. e# = —1), 

402 in which case T3 (resp. T2) describes a ferro-collincar 

403 configuration of the moments on each triangle. Inten- 457 

404 sity asymmetry of the pairs K ± f of magnetic satel- 458 

— * _ _ 45g 

405 lites about reciprocal nodes K indicated that for, a left- 

406 handed structural chirality €t = —1, if £h = — 1 then 

407 Ti is selected and if e H = +1 then T 2 is selected. This 

~* . . 462 

4os interdependence of the dephasing of moments within and 

. 463 

409 between the Bravais lattices C v was explained as arising 

410 from the twist in the exchange paths connecting the mo- 

411 ments of consecutive a,6-plancs, which depends on the 

412 structural chirality £t and imposes the magnetic trian- 

413 gular chirality ct^h- X-ray anomalous scattering con- 

414 firmed that the structural chirality of the investigated 

415 crystal is ct = — 1. Neutron spherical polarimetry finally 

416 demonstrated that only the magnetic helicity cr = —1, 

417 and therefore only the (£h,£t^h) = ( — magnetic 

418 helicity-triangular chirality pair, is selected, which was 

419 ascribed to the uniform Dzyaloshinsky-Moriya interac- 

420 tions with the Dzyaloshinsky-Moriya vectors all along the 

421 c-axis. This model was later confirmed by polarized neu- 

422 tron inelastic scattering with polarization analysis, which 

477 

423 allowed probing both the symmetric and antisymmetric 

424 nature of the dynamical correlations associated with the 

1 ^479 

425 magnon excitations emerging from the magnetic order. 

480 

426 A crucial point of the reported model of the circular 481 

427 helices with moments within the a,6-plane is that the 482 

428 dephasing of the moments associated with the triangu- 483 

429 lar configuration of moments on each triangle leads to 484 

430 zero magnetic structure factors at the scattering vec- 485 

431 tors (0,0, ±r). One however may recall that the neu- 

432 trons detect only the components of the moments per- 

433 pcndicular to the scattering vectors. An additional 

434 sine-wave amplitude modulated component along the c- 

435 axis of the moments is therefore not to be excluded, 

436 in which case we would rather have the combination^ 

437 v v (Tr/2,ct) v )e-^ + ae H v v (n/2,ct) v + vT/2)e-' i( «-- 7r/2) -U7 

438 v„(0, 0)e -i ^ = m a ^(x v +iotRy v )e~ % ^ v +m c z v e~ % ^<> . The4aa 

439 length m c of the vector v v (0,0) should however be smahV) 

440 enough so that the neutron intensities to which it should49o 

441 give rise at the other scattering vectors, (h, k, i ± r) withwi 

442 h 7^ or k 7^ 0, are drowned beneath the statistical un-492 

443 certainties of the neutron intensities associated with the493 

444 main helical modulation component. Table I actually il-494 

445 lustratcs that this c-component of the moments would495 

446 lead to a zero magnetic structure factor for the scatter-496 

447 ing vectors (0,0, ±r), and therefore would no longer be497 

448 detected by resonant x-ray scattering, if the stabilized ir-498 

449 reducible representation is either T2 or IV A non-zero499 

450 magnetic structure factor vectorially oriented along thesoo 

451 c-axis is computed only in the case of the irreducible rep-501 

452 resentation r x : F^ 1 = (0,0, f z ). The magnetic intensity^ 

453 I e ' e = F e , e F*, e (* stands for complex conjugation) in the503 



different diffraction channels 53 (e = er, ir and e' = a' ,tt'), 
associated with this amplitude modulated c-component, 
can be calculated with the help of Eq. (1) leading to 



= , 

l a ,K oc sin 2 6 B 



(2) 



where 9 B is the Bragg angle. Noteworthy is the absence 
of any azimuthal dependence. We therefore expect no 
modulation of the intensity as we rotate the sample about 
the scattering wave vector. Moreover, we expect l a = 
(la' a + I-rr'a) = Ijt = (I*'* + 1-x'tt) and Ir = I L . The 
latter equality can be derived from Eq. (Al) in Ref. 54 
which states I R - I L =lm{F^F a , a + F*, v F^ a }. 

Another deviation of the magnetic structure might 
arise from a slight ellipticity of the helices, but according 
to Table I this would remain invisible in the case of the ir- 
reducible representation T\. A finite magnetic structure 
factor, either F% = (f x ,f y ,0) or F% = (&,ft t 0), for 
the scattering vectors (0,0, ±r) would be obtained only 
if cither the T2 irreducible representation or the T3 irre- 
ducible representation were to be stabilized as the main 
helical modulation component of the magnetic structure, 
but this is ruled out from the neutron diffraction data. 

A mixing of the irreducible representation Fi with the 
irreducible representation T2 (or T3) finally is a priori not 
to be excluded, though this would imply that the mag- 
netic transition is necessarily first order. Nevertheless, 
the additional magnetic component should be extremely 
tiny to escape standard powder neutron detection, since 
it should lie in the a,6-plane to produce a non zero mag- 
netic structure factor. In the case of the ferro-collincar 
configuration in the a,6-plane, associated with irreducible 
representation T2 for en = — 1, which gives rise to a mag- 
netic structure factor of the form F^ 2 = (f x ,fy,Q), one 
calculates with the help of Eq. (1) the intensities: 



Ia'a = , 

I-rr'cr = la'TT = COS 2 I 

I^ = k 2 sin 2 (2fl B ), 



(3) 



where the constants k% depend on the amplitude of the 
component of the moments associated with the irre- 
ducible representation T2 and their orientation in the 
a,6-plane with respect to the moments associated with 
the main irreducible representation Y\. Even in this 
case there is no azimuthal angle dependence, but we find 
I a < In and Ir ^ II- Including both Ti and T3 con- 
tributions will lead to an azimuthal angle dependence in 
the rotated channels and again I a < 1^ and Ir =/= II- 

We are now in the position to compare the x-ray ex- 
perimental data with the prediction from representation 
theory. Fig. 2 shows that the ratio Ir over II is constant 
as a function of the azimuthal angle and equals one. Also 
the ratio I a over 1^ is roughly constant within the error 
bars and is very close to one. It is thus clear that no 
mixing of irreducible representations is detected and that 
the magnetic structure abides by only the irreducible rep- 
resentation Ti but involves components of the moments 
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504 along the three orthogonal direction in space. As a wholes54 

505 it consists of moments in a triangular arrangement onsss 

506 each triangle in the a,&-plane helically modulated alongsse 

507 the c— axis and exhibiting small up and down oscillations557 

508 along the c-axis in phase with each other and with thesss 

509 same period as the helical modulation. Such a motif iss59 

510 reminiscent of the beatings of butterfly wings (althoughseo 
5n these wings here are three in number and not four),56i 

512 that lead us to dub it as "helical-butterfly". The ex-562 

513 istence of the butterfly component is consistent with the563 
5u Dzyaloshinsky-Moriya interactions. Owing to the pres-564 

515 ence of the three 2-fold axes at 120° of each other in theses 

516 a,5-plane, each being perpendicular to one of the three 

517 sides of every triangle of moments, the Dzyaloshinsky- 

518 Moriya vector associated with each pair of moments must^ 
5w by symmetry lie within the plane containing the link con- 

520 necting the two moments. 55 The Dzyaloshinsky-Moriya 567 

521 vector field may therefore have a uniform component 568 

522 along the c-axis and a multi-axial component along the 569 

523 side of each triangle. It is this last component that gives 570 

524 rise to the butterfly component. It has been suggested 571 

525 that its contribution might be significant 56 if not domi- 572 

526 nating. 18 

527 Let us now analyze the azimuthal-angle dependence of 574 

528 the (0, 0, 2t) reflection. According to the Ti magnetic 5 75 

529 structure factor F^ 1 = (0, 0, f z ) and the formalism to cal-5 76 

530 culate magnetic diffraction intensity in Ref. 24 we should 57 7 

531 observe intensity only in the unrotated ir'ir scattering 578 

532 channel which is at odds with the data shown in Fig. 3. 579 

533 To reconcile the observations with theoretical prediction5 80 

534 we must adopt a more sophisticated model which does 58 i 

535 not rely on the fact that the resonant ion environment is 582 

536 cylindrically symmetrical. We need a tensorial structure 58 3 

537 factor &q where the positive integer K is the rank of the 58 4 

538 tensor, and the projection Q can take the (2K-\-\) integers 

539 values which satisfy —K<Q< K . For a dipole tran- 

540 sition, tensors up to rank 2 contribute (K < 2). K = 

541 reflects charge contribution, K = 1 time-odd dipole, and^ 

542 K — 2 time-even quadrupolc. For our superstructural re- 

543 flection we are interested in the quadrupolar contribution 587 

544 and given the presence of the 3-fold axis parallel to the 588 

545 c-axis we have *£(0,0,2r) = (-1) 2t (T£') [1 + 2 cos(Qa)] 589 

546 which is non-zero only for Q = 0. (Tq) is an atomic 

547 tensor that describes the contribution of each atom to 

548 the structure factor. Making use of the formula in ap- 590 

549 pendix C of Ref. 37 we obtain the following results for the 

550 structure factor in the different polarization channels: 

Ft = F > = 

Fir< n OC -=(1 + COS 2 <9b) *o > 

V6 

F a /F„ = -1/(1 + cos 2 B ). 

551 A derivation of such relations is presented in the Ap- 

552 pendix. Such a model suggests no azimuthal dependence 

553 in all the diffraction channels and a ratio I a j 'I v = 0.6 in 



relative agreement with the azimuthal dependence shown 
in Fig. 3 with a \ 2 = 6.1. Agreement is improved 
(x 2 = 2.2) by letting the ratio value vary as a free pa- 
rameter, with the experimental value of 0.54±0.02, still 
reasonably close to the one derived by Eq.(4). However, 
such a ratio, as exemplified in Fig. 1, is not constant as a 
function of energy. These deviations might arise from a 
small symmetry break resulting in a loss of the 3-fold axis 
which would cause extra terms to appear in the structure 
factor. Experimental uncertainties are however too big 
to extract more quantitative conclusions on the presence 
of such contributions. 



VI. CONCLUSION 

We have studied the magnetic structure of the in- 
triguing compound Ba3NbFe3Si20i4 with resonant x-ray 
diffraction at the Fe L edges and O K edge. These ex- 
periments give new insight into the details of the mag- 
netic structure recently determined by neutron diffrac- 
tion. Our experiments have found an extra sinusoidal 
modulation of the Fe magnetic moments along the crys- 
tallographic c-axis, concomitant with the helical order 
in the a,6-plane, generating an helical-butterfly magnetic 
structure. Such sinusoidal modulation arises from the 
Dzyaloshinsky-Moriya interaction as suggested by sym- 
metry consideration and recent linear spin-wave theory 
calculations. 56 The orbital magnetic signal observed at 
the oxygen K edge reflects the strong hybridization be- 
tween iron 3d and oxygen 2p states. Finally, the energy 
dependence of I a /In ratio for the (0, 0, 2t) reflection 
hints to a possible symmetry break with loss of the 3- 
fold axis, however ab initio calculation would be needed 
to obtain quantitative informations. 
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Appendix: Quadrupolar structure factor 



< = ^ (A.l) 



A\ = -±= - ttj) (A.2) 



591 In analogy with Ref. 37 we obtain expression for ^q, 

592 written in the coordinate space (x,y,z), as a sum of quan- 
(4)593 tities that are even (Aq) and odd (-Bq) functions of the 

594 projection Q with —K<Q<K. 

595 We give expression analog to Eq. (B5) of Ref. 37 for a 

596 generic (0,0,£) reflection: 
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598 Limiting ourselves to the quadrupolar contribution 

599 (K=2) and taking advantage of the structure factor 
/g x eoo *|(0,0,2r) = (-1) 2t (T^)[1 + 2cos(Qa)] we have only 

A 2 = — ($ 2 _ 2 + $1) - - V 2 (A.3)«n different from zer0 . 

602 Expressions in Eq. (A. 3) therefore simplify leading to 

A\ = - (^ 2 _ 2 - ^l) 603 e.g. B 2 Q =0 and A\ cx A\. Substituting Eq. (A.3) in 

Eq. (C1-C3) of Ref. 37 one obtains the expression quoted 



2 

B\ = i - $1) 



606 
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609 
610 
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613 
614 
615 
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